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Abstract 

The paper is devoted to the generalized wall functions of Robin-type and their 
application to near-wall turbulent flows. The wall functions are based on the trans- 
fer of a boundary condition from a wall to some intermediate boundary near the 
wall. The boundary conditions on the intermediate boundary are of Robin-type 
and represented in a differential form. The wall functions are formulated in an 
analytical easy-to-implement form, can take into account the source terms of the 
momentum equation, and do not include free parameters. The log-profile assump- 
tion is not used in this approach. A robust numerical algorithm is proposed for 
implementation of Robin-type wall functions to both finite-difference and finite- 
volume numerical schemes. The algorithm of implementation of the Robin-type 
wall functions to existing finite-volume codes is provided. The axisymmetric im- 
pinging jet problem is numerically investigated for different regimes on the base of 
the wall-functions implemented to the high-Reynolds-number k — e model. 

1 Introduction 

Problems related with turbulent near wall flows appear in many industrial applications. 
It is well-knoen that turbulence vanishes near a wall due to both the no-slip boundary 
condition for the velocity and the blocking effect caused by the wall. In the vicinity 
of the wall, there is a thin sublayer with predominantly molecular diffusion and viscous 
dissipation. The sublayer has a substantial influence upon the remaining part of the flow. 
An adequate resolution of a solution in the sublayer requires a very fine mesh because 
of the thinness of the sublayer and high gradients of the solution. It makes the model 
to be time consuming and often it is not suitable for real design. Because of the low 
turbulent Reynolds number in the sublayer, the models that resolve the sublayer are 
called low-Reynolds-number (LR) models. 

In turn, the high- Reynolds-number (HR) models do not resolve the viscous sublayer. It 
significantly saves computational efforts [1,2]. In the HR models, the boundary conditions 
or near-wall profiles are represented by wall functions. The wall functions usually are semi- 
empirical and have very limited applications [1-5]. First wall functions are based on the 
log-law profile assumption for the velocity [4,5]. In addition, their formal extension to 
complex flows demands time-consuming iterations for calculating the skin friction included 
in the log-law profile. A substantial disadvantage of these wall functions is a strong 
dependence on the near wall mesh used. This problem is especially pronounced if the 
first mesh point is located inside the viscous sublayer. To avoid this, the scalable wall 
functions are suggested in [6] . Wilcox assumes [7] that the pressure gradient must be taken 
into account to avoid the mesh dependence. Yet, the recently suggested adaptive wall 
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functions [8] overcome this limitation by using look-up tables for turbulent quantatives 
and skin friction. 

In more recent wall functions [1-3,9,10] source terms, such as the pressure gradient, 
might be taken into account. The numerical comparisons done in [1-3, 10] showed that 
such advanced wall functions give substantially better prediction than the standard wall 
functions. In [3,10], the analytical wall functions are obtained by approximate integrating 
boundary-layer-type equations in the wall vicinity using the assumption that all terms 
besides the diffusive one are constant. At the wall, the boundary conditions are the same 
as those used in the LR models. An analytical profile for the turbulent viscosity are then 
used in the cell nearest to the wall to reconstruct the near-wall solution. The wall functions 
for the turbulent kinetic energy and its dissipation are based on the local analytical 
solution for the velocity in the near-wall cell. In computations, the numerical flux to the 
wall is taken from the previous iteration as it is performed in the case of the standard 
wall functions. Although approaches [1-3, 9, 10] allow one to make substantially better 
prediction in comparison to the standard methods, their realization seems to be quite 
complicated. The wall functions [1-3,9,10] are only represented in a finite-difference form. 
Although this form is suitable for finite-volume algorithms, its use for finite-difference 
approximations is not clear. Similar to the standard wall functions, this approach faces 
substantial problems if the nearest to the wall cell is in the viscous sublayer. Also, it is 
important to note that the second to the wall cell cannot be much smaller or bigger than 
the first one because of the integration over the first cell. 

The method of boundary condition transfer is suggested in [13,14]. The method al- 
lows us to transfer a boundary condition from the wall to some intermediate surface. The 
boundary condition is transfered either approximately (analytically) or exactly (numeri- 
cally). The boundary conditions on the intermediate surface are always of Robin-type (or 
mixed type) and represented in a differential, mesh independent, form. These boundary 
conditions are set on both a function and its normal derivative. Therefore, their real- 
ization does not require additional iterations between, for example, the velocity and skin 
friction. This brings an additional robustness to the algorithm of their implementation. 
These boundary conditions are interpreted as generalized (Robin-type) wall functions. 
Another advantage of these wall functions is related with their universal formulation for 
all dependent variables. 

The Robin-type wall functions take into account the influence of the source terms in 
governing equations. The location of the point, to which the boundary conditions are 
transferred, does not make any considerable effect on the mesh distribution nearby this 
point. The wall functions can be implemented in both finite-difference and finite-volume 
approximations. It is shown how the wall functions can be implemented in existing codes. 
Preliminary tests for channel flow [15] and impinging jet [14] have shown promising results 
in terms of both accuracy and efficiency gains. In comparison to the analytical wall 
functions [3, 10], the key advantages of the Robin-type wall functions are related with 
their robust implementation and universal differential formulation. 

In the paper below the Robin-type wall functions are implemented in the k — e model 
and applied for the axisymmetric impinging jet problem. The computational results are 
compared against available experimental data. Implementation of the wall functions to 
both finite-volume and finite-difference schemes are discussed. 
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2 Generalized (Robin type) wall functions 



In order to formulate the Robin-type wall functions, first let us record the governing 
equations in the following general form: 

{HUy) y = R h (1) 

with Dirichlet boundary condition on the left-hand side: 

u (0) = n (2) 

Equation represents the general form of the boundary-layer-type equation. The 
right-hand side Rh is an appropriate source term including, e.g., the pressure gradient in 
the momentum equation. 

If the right-hand side Rh = const, then the Robin-type wall functions are formulated 
as follows [13, 14]: 

<v')=^h%(v')-^h, (3) 



where 







Relation ((HJ) can be interpreted as a boundary condition of Robin-type transferred 
from a wall (y = 0) to some point y*. This boundary condition can be either exact (if 
the exact function of \i is used in or approximate (if fi is estimated by one way or 
another). One should emphasize that the Robin-type boundary condition is set for both 
a function and its derivative. It is easy to see that the mesh distribution nearby the 
point y* can be independently chosen on the location of this point. Implementation of 
Robin-type conditions to both finite-difference and finite-volume schemes are considered 
below. 

In the general case Rh = Rh(y), we have 



where 







Having assumed that the coefficient varies piece-wise linearly 



if < y < y v 



it is possible to obtain analytical expressions for fi and /2 if Rh = const and y v < y* : 

fi = a IJ y v {l + 9hia lt ), f 2 = a„y v [(l - 6)y* + y v {6 2 a^ lna M - 1/2 + 6)} , (7) 
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where a u = a*/a w ,6 1 = ^ — ^ . The parameter 6 1 represents cotangent of the 
inclination angle of the dependence \ij fi w on y/y v . 

If Rh = const, the wall-flux can be found as follows [14]: 

r w = ^f^- + {h/h - y*) R h . (8) 
h 

This formula can be easily generalized on the case of a variable right-hand side Rh 
using (jSJ). 

The method of boundary condition transfer technique can be used to derive the wall 
functions for the tangential and normal velocity components U and V, temperature T , 
and turbulent kinetic energy k. 

Having neglected diffusion parallel to the wall, the momentum and enthalpy transport 
equations can be written in the Cartesian coordinate system (x, y) as follows: 



d_ 

dy 
d_ 

dy 



,dU 
dV 



d_ 

dy 



Pr 



dT 
Pr t dy 



pU— + p V— + — 

dx dy dx 

TT dV jr dV dP 
pUj- + pV— + — 
dx dy dy 

pU dT +p ydT 

dx dy 



(9) 
(10) 

(11) 



Here fii and fi t are the laminar and turbulent viscosities, accordingly; Pr and Pr t are 
Prandtl numbers; U and V are the velocity component in the (x, y) coordinate system; p 
is the density; P is the pressure. 

The intermediate boundary conditions for U, V and T at point y* are given by (0) 
upon substitution U, V or T instead of u accordingly Evidently, the coefficient ji in (JT]) 
must be considered as either fii + p t or fii/Pr + fi t /Pr t . In the case of the momentum 
equation uq = 0. If y* is chosen in the vicinity of the wall, the right-hand side Rh can be 
simply evaluated at y* . Thus, in the case of the momentum equations and enthalpy the 
relative right-hand sides follows: 



R h = R hu = P (u-( y *) + V-( y *)) + -(y* 
( dV dV \ dP 

dT dT 
Rh = Rht=p(U—(y*)+V—(y*)) 



(12) 
(13) 
(14) 



Thus, all the terms of the parabolized (reduced) Navier-Stokes equations (PNS) [16] are 
taken into account. It worth noting that in the boundary condition © for the normal 
velocity it is not assumed to be zero. Thus, these wall functions do not have direct 
restrictions to their exploration in modeling separated flows. 
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Unlike [3], a similar approach is applied to the equation for the turbulence kinetic 
energy as well: 



d_ 

dy 



fit sdk 



where Pk is the production of the turbulent kinetic energy, e is its dissipation; Pr^ is the 
Prandtl number. 

Having evaluated the convective terms, we obtain the following expression for the 
right-hand side R^. 

( dk , dk . / dU\ 2 

R k = R hk = p (U-( V -) + V-W) ) + pe - p t (-J (16) 

Having assuming a piece- wise linear behavior of the function fi t '- 

Jo, H0<y<y v , . 

* iiy v <y<y*, {n} 

where y v is the thickness of the viscous sublayer near the wall, the coefficients fi and fi 
in (jHJ (the latter term only if Rh = const) can be evaluated by (jJJ). 
For the momentum equation 

<*n = H / Pi, = -, /x = A*t + A*ti ( 18 ) 

while in the case of the energy equation 

fii/Pr + fi* t /Pr t y* - y v Pr t m 

a v = T5 1 & = 5 * ( 19 ) 

If the turbulent viscosity /x* in (fTTj) is evaluated as follows [3]: 

/i t = C^Cip y v hi = C^CiRe v /i; « 2.5 // z , (20) 



where fc* = k(y*), C = 0.09, Q = 2.55, = = 10.8, then is a constant equaled 

to 0.4 for the momentum equation. 

It has been found from the computations that it is more accurate to evaluate the 
turbulent viscosity at y* from the HR k — e model directly 

Ht = C,p(k*)) 2 /e (21) 

rather than from equation (|2*UI) . 

The sublayer thickness y v is evaluated as follows [3]: 

y v = Re v fn/ (p-\/k] , (22) 

where k v is the value of k at the edge of the viscous sublayer. 
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If y* < y v , then the boundary conditions are set inside the sublayer, and formulas J7J) 
are not formally valid. It is suggested to pose the boundary conditions at the edge of the 
sublayer y = y v as in [6] because y v is small enough. Then, the coefficients f\ and f 2 in 

can be evaluated as follows: 

fi = a^yv, f2 = oi^yl/2. (23) 

It is then assumed that the turbulent viscosity fi t reaches value ()21|) at the edge of the 
viscous sublayer immediately. These boundary conditions are consistent with boundary 
conditions (J7J) taking in the limit y* — > y v or 9 — > 0. 

The dissipation of the turbulent kinetic energy e is evaluated as in [3]: 

I Ciy ' else > 

where y d = 2Cifn/ (j?Vk*j ■ 

The wall function for the turbulent energy k is used in form Q, © and depends 
on the right-hand side Rhkiy) represented by equality (JEJ). It includes the dissipation e 
and derivative dU/dy. The former term is taken from (|2*^ while the latter term can be 
evaluated in the interval [0,y*] from the reduced momentum equation ©, (fT2*|) : 

+ iM)dU/dy = [(fit + fi)U{y*) + f 2 R hu ] / h + (y ~ y*)R hu , (25) 

where the turbulent viscosity /i t is defined by fITTj) . Equation (|23j) is obtained by the 
integration of equation Q and use relation (J3J) to exclude dU/dy(y*). Thus, this equation 
relies on the PNS assumptions used. 

To evaluate y v from (J22J), it is possible to use the value k* instead of k v . It allows us 
to simplify the evaluation algorithm for y v substantially. First, a similar opportunity was 
noticed in [3]. The estimation of Re v is varied between 10.8 and 20 [3,18]. It corresponds 
to the interval between the upper limit of the viscous sublayer and the point at which the 
linear and logarithmic parts of the velocity profile intersect for the channel flow [19]. It 
is not clear which value in this interval is most appropriate to approximation (fTTj) . In all 
computational results given below, Re v = 12. 

It worth noting that the coefficients j\ and f 2 in wall functions (jUlH© depend only 
on y* and k*. The latter value is determined from the solution of the HR model at the 
boundary point y*. Hence, the intermediate boundary conditions (jSJ) at y = y* complete 
the boundary- value problem in the interval [y*,y e ] and can be considered as generalized 
wall functions. These boundary conditions are of Robin-type and similar to the "slip 
boundary condition" at the edge of the Knudsen-layer in aerodynamics. One should note 
that the boundary conditions are linear with respect to the leading variable. As it follows 
from Q and (0), the source terms in the wall functions can only be essential far enough 
from the wall because of the quadratic dependence of f 2 on y*. 

The HR solution obtained in [y*, y e ] can be extended to interval [y v , y*] using the 
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analytical solution in this interval: 



u(y) 



u(o) + MvK - My) 

a„y v (l + 9\nn(y)), 

&0V [y* -9y + (^ 2 « M " 




(26) 





0: 



'2 



1 



1/2 + 9)y v ] , 




Thus, the intermediate boundary is not necessarily to be related to the nearest to the 
wall cell. It is possible to take y* far enough from the wall and complement the solution 
on the region of the sublayer by (f2*oT) . 

It worth noting that, although the problem is solved in the bulk domain [y*, y e ] only, 
the flux to the wall (e.g., skin friction) can be easily evaluated considering (JHJ) (or its 
analogy for the temperature in the case of heat flux). 

Thus, the developed wall functions can be applied to all dependent variables but e in 
a uniform manner. It is possible to show that the coefficients f\ and f2 can be determined 
analytically even in the case of the turbulent kinetic energy k. This fact can be useful 
for saving computer resources. For engineering purposes it worth considering a simplified 
version of the wall functions corresponding to j'2 = 0. Since the coefficient fi can be 
chosen to be the same for all variables, the wall functions become fully identical in this 
case. 

The generalized Robin-type wall functions are not based on a numerical approxima- 
tion in the inner region [0,y*], as in [1-3,10], therefore the location of the intermediate 
boundary is not very substantial for the mesh distribution in the bulk domain. It means 
we can choose, e.g., a fine mesh despite a relatively big value of y* (or vice versa) without 
loose of stability. 

In this paper we make the main focus on the analytical evaluation of the coefficients 
fi and f2 in (J3J). Yet, the coefficients can be evaluated numerically by integrating LR 
equations in the interval 0, y*. It may lead to the decomposition method [13,14]. Also, 
this approach can be naturally integrated in the numerical wall-functions [2,9]. 

3 Numerical implementation of Robin— type wall functions 

The Robin-type wall functions can be implemented to both finite-difference and finite- 
volume RANS approximations. In this section, some aspects of their robust implementa- 
tion are considered. 

A general remark related with robust implementation is as follows. Boundary condi- 
tions of Robin-type are set on both a function and its derivative. Upon approximating the 
derivative, both terms should be considered at the same iteration (or time step). Taking 
into account one of the terms from a previous iteration leads to additional iterations, at 
least. It is easy to see this property in the case of a linear equation. Only simultaneous 
consideration of both terms provides an iterationless solution. 

The boundary condition is represented by (JHJ) for any y* in the vicinity of the wall 
even if y* vanishes. This boundary condition fully replaces the original boundary condition 
u(0) = u . 
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In numerical simulation of turbulence, the finite-difference numerical schemes pre- 
serving positiveness of a solution [17] are very efficient because unknown variables such 
as the turbulent kinetic energy k or its dissipation e must be positive. The following 
numerical procedure can be used for developing the positive definite schemes in solving 
boundary- value problems with Robin-type boundary conditions [14]. 

Boundary condition (jSJ) can be rewritten in the following general form: 

fc(0) = adk/dy(0) + (3, (27) 

assuming that both the function k and its derivative dk/dy are positive. This assumption 
is valid in the case of real physical problems for the turbulent kinetic energy in the wall 
vicinity. The coefficient a is positive because fi is always positive but the coefficient (3 
can be negative (mostly, where e > P&). In computations it can lead to a negative value 
of k. To avoid such a case, it is suggested to rewrite (J2*7)) in the following form if /3 < 0: 

fc(0) = adk/dy{0) + (3- k 



k-(oy 

or 

fc(0) = adk/dy(0), (28) 

where a = i_pJ L k -^ and k~(0) is the value of k(0) taken from the previous either time 
step or iteration. 

At last one should note that at some distance from the wall the derivative dk/dy 
becomes negative. Yet, in this area the function j3 is positive (P& > e) and, therefore, 
there are no principal difficulties in reaching a positive solution. 

Though Robin-type boundary condition (jHJ) can be set at the wall, its implementation 
to existing codes based on finite-volume schemes is more easy in the following treatment. 
Assume that unknown variables are defined at the centers of cells. For the sake of sim- 
plicity let us consider ID approximation in the normal to the wall direction. It is enough 
to consider approximation at the nearest to the wall cell since the rest approximation is 
remained without any modification. Let us denote values at the centre of the cell by index 
1/2 and values at the edge, opposite to the wall, by index 1. Then, considering y* = yy 2 
we have a relation in the following form: 

du ~ 

u 1 / 2 = u w + f 1 — + f 2 , (29) 

dy\i/2 

/•»l/2 R d 

where f 2 = -i^™/ 2 . At the first cell the governing equation JI} is then integrated 
only from yy2 to y\\ 

du du f Vl 
\i— - n— iti/2 = / R h dy (30) 
dy\i dy {1/2 J yi/2 

Equation (|3U|) represents the approximation of the governing equation in the near-wall 
cell. Here, the flux is approximated by an ordinary approach while the flux at yi/ 2 

is obtained from (|29jl . A typical approximation can be written as follows: 

U 3/2 ~ Ml/2 U 1/2 ~U w -f 2 f yl 

Aii m/2 7 = / Rhdy (31) 

Vl H Jy 1/2 
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This kind of approximation is used for all variables, but e, including k. In the last case, 
the right-hand side is rapidly changed and the integral in the right-hand side must be 
evaluated accurately enough. It can be easily done either numerically or analytically 
using the analytical expression for the velocity gradient (|25jl . 

If y* > y v , the considered above approximation is close to the approximation based on 
the analytical wall functions (AWF) [3]. The principal difference is related with robust- 
ness. In the case of the Robin-type wall functions both the function and its derivative 
(flux) are simultaneously taken into account while in realization of the AWF the flux is 
calculated in iterations. This difference becomes more clear in the case of a linear equa- 
tion. Then, the Robin-type formulation of the boundary conditions does not require any 
iterations. 

At last, a brief comment can be done with regard to a staggered mesh. In this approach 
the velocity is defined at the vertexes of a cell. In our consideration this means u\ is known 
instead of iti/2- In this case, the simplest way to remain the uniform approach is based 
on the Taylor expansion: 

du 2/1/2 d 2 u 
ui=u 1/2 + y 1/2 — +^--r~ 2 • (32) 
"2/ 1 1/2 z a V 11/2 

This yields a required relation between U1/2 and u\ since the first derivative takes place 
in (j2HJ) and the second derivative can be evaluated via the right-hand side Rh- Thus, the 
flux at the nearest to the wall cell is as follows: 

du Ux-u w h 2 /8 - f 2 

d v\i/2 h + h / 2 h + V 2 

The same technique with slight modifications can be used for implementation to un- 
structured codes. 



4 Impinging jet 

The problem of impinging jet appears in many industrial applications related with ei- 
ther heating or cooling processes. The heat transfer problem on a turbulent circular jet 
impinging onto a flat is well studied experimentally [11,20-22]. Therefore, it became a 
widespread test case for different turbulent models including the LR [12, 24-28] and HR 
k - e models [6,9,24,29]. 

The problem is formulated as follows. A fully turbulent air jet, generated in a pipelike 
nozzle, impinges on a flat surface at the right angle. The Reynolds number, based on 
the nozzle diameter D and the bulk velocity, equals to Re = 23000 and Re = 70000. 
The distance L between the nozzle and the surface is varied between 2D and 14.D. Air 
is treated as an ideal gas and considered under normal conditions at temperature 293A'. 
The heated surface has constant temperature T w = 314. 9K. 

The computational domain spans 13D in the radial direction. The grid includes 
150x100 (axial x radial) nodes and 150x200 nodes. For the validation purposes, pre- 
liminary comparisons of the results obtained on different meshes were done to check grid 
sensitivity. The boundary conditions at the edge of the nozzle are specified using the 
profiles for a fully developed turbulent pipe flow. The computations of the local Nusselt 
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number are done for the different values of y* or Re y * = pyk*y* / fii calculated at the 
stagnation point. 

Linear eddy-viscosity model (EVM) drastically overpredict the turbulent kinetic en- 
ergy in the stagnation point region by an order of magnitude [29] . It inevitably leads to 
the considerable overestimation of the heat flux. As a result, the linear LR k — e models 
give unacceptable overprediction by a factor of two, even more [2,26-28]. Furthermore, 
the linear k — e model, as well as other EVM, is not entirely justified around the stagna- 
tion point because of the anisotropy of the flow. To improve prediction, along with the 
non-linear EVM, some modifications of the EVM are used including the implementation 
of a realizability constraint [28], introduction of empirical formulas for the Prandtl num- 
ber [26] and heat flux [6]. The application of the wall functions, in most cases, is also 
showed a poor performance [2,24,29]. More or less reasonable prediction was achieved 
in [2] using the Chieng - Launder wall function [30] and the scalable wall functions [6]. It 
is to be noted that in the latter case the empirical correlation was used for the local heat 
flux. The generalized wall functions were applied in [14] to simulation of the impinging 
jet at Re = 23000 with L/D = 2 and L/D = 6. Apart from the heat flux, the prediction 
of the wall friction is considered in [14]. As was noted above, the overprediction of the 
heat flux was obtained in the computational solution. 

The same effect but more expressive is observed at Re = 70000. In FigQ the com- 
putational results are compared against the experimental data for L/D = 2. Here and 
below, the local Nusselt number is scaled by Re°' 7 Pr 0A where Pr = 0.9. The solution [12] 
based on the low-Reynolds number k — e model predicts a substantially higher heat flux 
than the high-Reynolds number model. This effect was obtained in many other publi- 
cations including cited above. The solution based on the wall functions corresponds to 
Re y * = 109. It is to be noted that the dependence of the solution on the parameter y* 
is quite weak. In the next example with L/D = 4, the solutions corresponding to differ- 
ent values of Re y * are shown in FigEl Though the value of y* is varied by an order of 
magnitude, the curves are quite close each other. 

The mean velocity profiles divided by the bulk velocity are shown in Fig|B]for L/D = 6. 
The experimental data are represented by square symbols while the computational results 
are shown by the curves. At the region of the low mean velocity nearby the axis of 
symmetry (r/D = 0.5) the prediction of the velocity is quite accurate. At r/D = 3, 
where the flow is decelerated, the prediction is not so good. At this location, substantial 
underprediction of the velocity in the near wall region and overprediction in the outer 
region were earlier noted for both the LR and HR linear k — e models [6, 12,27]. 

In Fig|3]the distribution of the local Nusselt number is shown for Re = 23000. It is 
given a comparison between the computational results and experimental data for L = 10D 
and L = 14D. In these examples the wall is located far enough from the nozzle, and in 
the computations the overprediction of the heat flux nearby the axis of symmetry is not 
observed. 

5 Conclusion 

The wall functions are formulated as boundary conditions of Robin-type and represented 
in a differential form. These wall functions take into account source terms. The wall 
functions are obtained in a compact easy-to-implement analytical form and they do not 
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include any adjustable parameters. The mesh distribution inside the computational do- 
main can be chosen independently on the location of the intermediate boundary. The 
implementation of the wall functions is robust due to their simultaneous formulation for 
both a functions and its normal derivative. The Robin-type wall functions are written in a 
universal formulation applicable to all dependent variables but e including the kinetic tur- 
bulent energy and normal velocity. General approaches to implementing the Robin-type 
wall functions to finite-volume and finite-difference approximations are suggested. 

On the base of the generalized wall functions implemented to the k — e model the 
axisymmetrical impinging jet is investigated. The computational results show a reasonable 
correspondence to the experimental data and weak dependence of the solution on the 
distance from a wall where the boundary conditions are set. 

Further research can be devoted to application of the wall functions to separated flows 
and extension of them to LES. 
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Figure 1: Local Nusselt number for the impinging jet. Comparison between HR, LR 
solutions and experiment for Re = 70000 and L/D = 2. 




Figure 2: Local Nusselt number for the impinging jet. Comparison of computational 
solution for different y* against experiment for Re = 70000 and L/D = 4. 
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Figure 3: Mean velocity in the impinging jet at r/D — 0.5; 3. Comparison of computa- 
tional solution against experimental data for Re = 70000 and L/D = 6. 




Figure 4: Local Nusselt number for the impinging jet. Comparison computational solu- 
tions against experimental data for Re = 23000 and L/D = 10; 14. 
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